Search for the Decay Kl — > iftvv using n° — > e + e 7 
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We report on a search for the decay Kl — > ir°vv, carried out as a part of E799-II, a rare Kl 
decay experiment at Fermilab. Within the Standard Model, the Kl — > tPvv decay is dominated by 
direct CP violating processes, and thus an observation of the decay implies confirmation of direct 
CP violation. Due to theoretically clean calculations, a measurement of B(Kl — > 7r c Vi>) is one of 
the best ways to determine the CKM parameter r\. No events were observed, and we set an upper 
limit B(K L -> tPvv) < 5.9 x 10~ 7 at the 90% confidence level. 
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The decay Kl — > ir°vv is dominated by direct CP vio- 
lating processes within the Standard Model through sec- 
ond order diagrams such as Z penguins H. Indirect CP 
violating and CP conserving contributions are expected 
to be highly suppressed ^-|| for the following reasons. 
First order decay diagrams, which lead to relatively large 
indirect CP violation in Kl — * 7T7t, do not contribute to 
Kl s — ► t^vv because of the absence of tree level flavor 
changing neutral current. Indirect CP violating contri- 
bution via second order diagrams is suppressed by five 
order of magnitude (e 2 ). Long-distance indirect CP vio- 
lating and CP conserving contributions from Kl — > 7r°7* 
and Kl — > 7r°7*7* intermediate states, which are signif- 
icant in Kl — y 7r°e + e _ and Kl — > n fj, + fx~, do not exist 
in Kl — * because the neutrinos in the final state do 
not couple to virtual photons . 

Following the Wolfcnstcin paramctrization of the 
CKM matrix B(K L — > iftvv) is proportional to 

i] 2 . The uncertainty of the hadronic matrix element in 
Kl ifivv is eliminated by the experimental measure- 
ment of T(K + — > 7r°e + ^) and the lifetime of Kl, which 
leads to an uncertainty of ±1.5% in the expectation of 
B(Kl — » tt°vP). In addition, due to the small uncer- 
tainty (~ 3%) in the next-to-leading order QCD correc- 
tion P|, B(Kl — > ifivv) gives direct access to r\. The 
current knowledge of the CKM parameters [pi allows us 
to predict B{K L 7r°w) to be (1 - 5) x 10^ 0. The 
uncertainty comes directly from the input CKM param- 
eters. As the theoretical calculations are unambiguous, 
an observation of the decay Kl — > ifivv at the sensitiv- 
ity of ~ 10~ n would indicate the existence of direct CP 
violation, and an observation outside the predicted range 
would indicate new physics jO]. 

It is experimentally difficult to search for Kl — > t^vv 
because the signature is only an isolated 7r°. The cur- 
rent upper limit, B(Kl — > n Q vv) < 1.6 x 10~ 6 at the 
90% confidence level, was obtained by using ir° — * 77 de- 
cay [[l3]]. We report on the search for Kl — > ifivv in the 
Dalitz decay mode (ir° -> e+e~7, 7r£>) with the E799-II 
experiment using the KTeV detector at Fermilab. The 
data were collected in 44 days of running in 1997. Using 
the Dalitz decay of 7r°'s made it possible to reconstruct 
decay vertex position of the 7r°'s, allowing us to measure 
Pt, the 7r°'s momentum transverse to the Kl beam di- 
rection. The p t played an important role in background 
suppression. 

Figure |l| shows a plan view of the KTeV detector. The 
elements of the detector relevant to this search are de- 
scribed below. Kaons were produced by 800 GeV proton 
beam, with a typical intensity of (3.5 ~ 5.0) x 10 12 pro- 
tons per 19 sec beam pulse, that struck a 30 cm long 
beryllium-oxide target at a targeting angle of 4.8 mrad. 
In the first (second) half of the run period, two neu- 
tral side-by-side beams with a solid angle of 0.25 /zsr 
(0.35 fisr) each were defined by collimators downstream 
of the target. A 7.6 cm long lead absorber was placed 



to reduce photons in the beams. A series of sweeping 
magnets removed charged particles in the beams. The 
two beams entered a 69 m long evacuated decay volume 
starting 90 m from the target. The vacuum was kept at 
10~ 5 ~ 10~ 6 torr. The downstream end of the volume 
was sealed by a vacuum window made of Kevlar and My- 
lar. The thickness of the vacuum window assembly was 
0.0035 radiation lengths (X ) in total Q. The neutral 
beam was mainly composed of neutrons and Kl$> with 
other long lived neutral particles, such as A°'s and 5 s. 
The relative ratios of neutron, A , and S° to Kl at the 
beginning of the vacuum decay region were measured to 
be 3.5, 0.02 and 7.5 x 10 -4 , respectively. The average 
kaon momentum was 70 GeV/c. Approximately 3% of 
the kaons decayed inside the vacuum decay region. 

The position and momentum of charged particles were 
measured using a spectrometer consisting of four drift 
chambers, two upstream and two downstream of a dipole 
analyzing magnet. The magnet had a momentum kick 
of 205 MeV/c. Each chamber consisted of two orthog- 
onal views (x and y), and had approximately 100 [im 
single-hit position resolution per view. An electromag- 
netic calorimeter with dimensions ofl.9m x 1.9m and 
27 X in depth was used for photon detection and elec- 
tron identification [jl5| . It was composed of 3100 pure Csl 
crystals. The calorimeter had two 15 x 15 cm 2 holes lo- 
cated near the center of the array to allow neutral beams 
to pass through. The energy resolution of the calorimeter 
was below 1% averaged over the electron energy range 2 
to 60 GeV. A scintillator hodoscope was placed just up- 
stream of the calorimeter for charged particle triggering. 
There were 8 Transition Radiation Detectors (TRD's) 
between the spectrometer and the trigger hodoscope for 
e/n separation. The TRD's consisted of polypropylene 
fiber mats as radiators and active MWPC volumes. 

The hermetic photon veto system detected photons 
missing the fiducial area of the calorimeter. The system 
consisted of 3 sets of counters: perimeter vetoes (PV) 1- 
9, a collar veto (CV) and a beam hole veto (BHV). Each 
photon veto counter had a sandwich structure of lead 
(tungsten in the CV) and scintillator. The total depth of 
radiator was 16 X Q for PVs, 8.6 X for CV, and 30 X a 
(equivalent to ~1 nuclear interaction length) for BHV. 
The PV covered the outer part of the calorimeter and 
the fiducial volume. The CV was placed just upstream 
of the calorimeter and around the two beam holes. The 
BHV was located downstream of the calorimeter and in 
the neutral beam region. The BHV was segmented into 
two transverse sections (one per beam) and three longitu- 
dinal sections (10 Xq each). The first longitudinal section 
was designed to detect photons, and the last one to de- 
tect neutrons. Downstream of the calorimeter, there was 
a 10 cm lead wall followed by a scintillator plane (hadron 
veto) to reject charged pions. 

The trigger was designed to accept events with two 
electrons and a photon so that Kl — * it^vv and 
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Kl — > e + e _ 7 decays were accepted. The Kl — ► e + e~7 
decays were used to measure the number of decayed 
Kl's. The trigger hodoscope and drift chambers were 
used to select two charged track events. The calorimeter 
was required to have an energy deposit greater than 18 
(24) GeV in the first (second) part of the running period. 
Events with significant energy in the photon or hadron 
vetoes were rejected. Events with three or four clusters 
in the calorimeter with a minimum energy of 1 GeV were 
selected by the hardware cluster counting system Jl6| . 
The TRD pulse height information was used to identify 
electrons at trigger level. 

The strategy in offline selection was to identify ir^ de- 
cays by reconstructing the invariant mass (m ee7 ) and se- 
lecting high p t events in order to suppress backgrounds. 
The pt cut was used because 7r°'s from Kl — > tt°vV have 
a higher kinematic pt limit than those from most of back- 
ground processes. 

In order to avoid human bias in the determination of 
selection criteria, a blind analysis was performed. A 
masked region was defined in the pt vs m eei plane as 
125 < m ee7 (MeV/c 2 ) < 145 and 160 < p t (MeV/c) < 
240. Monte Carlo (MC) simulation was used to optimize 
all cuts while data within the masked region were hidden. 

The offline event selection began with the identification 
of Wjj decays by requiring 125 < m ee7 (MeV/c 2 ) < 145 (~ 
±3ct). There were five categories in the remaining back- 
grounds such as: Kl — ► n ± e T i' (K e s); Kl — > ir + ir~ ir^ ; 
hyperon decays; Kl — > tt ^ and Kl — * 7r 7r 7i£, ; beam 
background. Below we describe the cuts to suppress each 
background. 

One serious background was K e z decays where a pho- 
ton was radiated from the electron or overlapped acci- 
dentally, and the pion was misidentified as an electron. 
Electrons were selected by requiring 0.95 < E/p < 1.05 
where E is the energy deposited in the calorimeter and 
p is the momentum measured by the spectrometer. This 
cut was 94% efficient for detecting both electrons and 
0.4% for a pion. The transverse shower shape at the 
calorimeter was also used to distinguish electrons from 
pions. The confidence level to identify pions formed from 
the 8 TRD's was required to be less than 1%, which gave 
a 95.0% efficiency for electrons. Events with out-of-time 
accidental energy in the calorimeter were rejected. The 
photon energy was required to be greater than 3 GeV 
because accidental and radiated photons typically have 
lower energy. Dalitz decays, which favor low m ee , were 
selected by requiring m ee /m ee7 < 0.3, where m ee is the 
invariant mass of the electron pair. Defining 9+{6-) as 
the angle between a photon and a positron (electron) in 
the kaon rest frame, cos#+ + cos#_ was required to be 
less than -1.5, because tt and in semileptonic decays 
prefer a wide opening angle, thus a peak of zero in the 
cos + cos 9- distribution. These two kinematic cuts 
rejected 99.6% of K e ^ events with a signal efficiency of 
78%. 



Backgrounds involving ir^ decays with unrecon- 
structed charged particles, such as Kl 7r + 7r~7T£, , were 
suppressed by eliminating events with more activity in 
the drift chambers than expected from two charged track 
events. 

High momentum A°'s and S°'s could reach the decay 
region in spite of their short life time. Decays of these 
hyperons could lead to backgrounds such as A — » mr Q D 
and S° — > A ^, because of undetected neutrons. These 
backgrounds were reduced by requiring the z position, or 
decay distance from the target, to be greater than 120 m. 
Since hyperons decaying in the decay region had higher 
energy than kaons, typically 200 to 300 GeV/c, events 
with photon energy greater than 50 GeV were rejected. 
In order to suppress backgrounds with neutrons such as 
A — > nn° D , the energy deposited to the third segment of 
BHV was required to be less than 200 minimum ionizing 
particles equivalent. This cut was applied only for the 
+x(— x) side of BHV when the decay vertex was found 
in the +x(~ x) region to minimize the signal loss due to 
accidental activity. 

The Kl — > ^ir^ and Kl —> ir ^ ^ backgrounds 
were suppressed by the photon veto system. The thresh- 
olds for measured photon energy were set to 200 MeV for 
PV1 and PV2, 250 MeV for PV3, 100 MeV for the rest 
of PV's, 1 GeV for CV, and 5 GeV (8.5 GeV) for the first 
section of BHV on the same (opposite) side as the recon- 
structed decay position. The electromagnetic calorimeter 
was also used as a part of the photon veto system. The 
number of clusters with energy greater than 1 GeV was 
required to be three, and events with extra clusters with 
energy greater than 250 MeV were rejected. The pho- 
ton veto requirements rejected 99.5% of Kl — ► and 
over 99.99% of K L -> 7r 7r 7r°, events, while 41% of the 
signals, as measured by Kl — ► e + e _ 7 events, was lost. 
Most of the signal loss arose from the BHV, because the 
high rate neutral beams (13 MHz Kl and 44 MHz neu- 
tron) struck the BHV. 

Another background was associated with 7T°'s pro- 
duced by beam interactions with detector materials, pri- 
marily the vacuum window. Figure |^ shows the pt vs 
z distribution for data. There was a cluster of events at 
z ~ 159 m, the location of the vacuum window. To reject 
such events, the decay vertex position in z was required 
to be less than 150 m. 

The remaining backgrounds were primarily from hy- 
peron decays, which had a well-reconstructed ir^ decay 
in the fiducial region. These were rejected by requir- 
ing p t to be 160 < pt(MeV/c) < 240 as shown in Fig- 
ure |^. The cut on the high end was determined from 
the kinematic limit of Kl — * ifivv decays, allowing for 
resolution. The main peak arose from A — > nit® D , and 
the shoulder at 135 MeV/c was from S° -> A vr^. The 
MC events were normalized by the absolute number of 
decayed Kl's, A°'s and S°'s. With this absolute nor- 
malization, the agreement between data and MC distri- 



3 



butions is excellent. Combining m ee7 and p t cuts, the 
efficiency for A — > mr^ was less than 1.4 x 10 -6 , and 
O(10 -2 ) to O(10~ 5 ) for the other S° decay backgrounds, 
while the signal efficiency was 46%. 

In order to verify the MC simulation and our under- 
standing of the backgrounds, events around the masked 
region were compared between data and MC as shown in 
Figure [|. The region (f), which had the largest discrep- 
ancy of all the regions, had a Poisson probability of 5.6% 
for observing 10 events when 6.5 events were expected. 
The good agreement between the expectation and the 
data in both the pt shape and the number of events vali- 
dates the Monte Carlo simulation and our understanding 
of the backgrounds. Even if one or more of the cuts is 
relaxed, the agreement is still excellent. 

Using the MC, the background levels except for those 
from beam interactions were estimated and summarized 
in Table Q. In the case of backgrounds associated with 
beam interactions, the sideband data were used for the 
estimation. Using the shape of the tail in the z distribu- 
tion for the p t > 240 MeV/c region, as shown in Figure ||, 
the contamination of backgrounds to the signal region 
was expected to be 0.04 events (no correlation between p t 
and z was found). In total, 0.12±°;^ background events 
were expected. 

The signal acceptance for K^s decaying between 90 m 
and 160 m from the target and with a momentum range 
of 20 to 220 GeV/c was calculated from MC to be 0.152%. 
The acceptance for Kl — > e + e~7 was similarly calculated 
to be 0.815%. The single event sensitivity (SES) of this 
search can be expressed as: 
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where A, N, and B represent the acceptance, the num- 
ber of observed events and the branching ratio of each 
mode, respectively. Based on the above equation with an 
observed number of Kl — > e + e _ 7 events of 15951 and a 
branching ratio of 9.1 x 10~ 6 |fL7|] , the SES was calculated 
to be [2.56 ± 0.02(stat.) ± 0.17 (ay a.)] x 10" 7 , where the 
sources and sizes of the errors are summarized in Table ||. 
The uncertainties on branching ratios of Kl — » e + e~7 
and 7r^, contributed a large part of the systematic error. 

Finally, we examined the signal region and found no 
events. Since no signal events were observed, the 90% 
confidence level upper limit on the branching ratio was 
determined to be B(K L — > ir a vv) < 5.9 x 10~ 7 . 

We dehne "background limit" as SES multiplied by the 
number of expected backgrounds. This figure of merit 
shows the experimental potential for rare decay searches 
because it takes account not only a SES but also an ex- 
pected background level. The "background limit" in this 
search is 3.1 x 10 -8 . This is a factor of 49 lower than the 



"background limit" of reference |ll|, in which the SES 
is 4.04 x 10~ 7 with an expectation of 3.7 background 
events. Even with this better technique, direct CP viola- 
tion has not been observed in this mode. To observe the 
Kl — > 7r°z/P at the Standard Model level, a search with 
four orders of magnitude better "background limit" will 
be required. 
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FIG. 1. The KTeV detector configuration for E799-II. 
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TABLE I. Summary of expected background contribution 
in the final signal region. 



Decay mode 


Expected number of events 


Kl —> nev ± 7 




0.02 ± 0.02 


K L -> n + TV'TT% 




< 0.01 


A — ► n-K° D 




< 0.04 


H°-» A (-» JOT- 
S' 1 -f A°H Tin 
S° -» A°H nvr^ 


)A 
A 


n m +0.006 
u ' ul -0.004 
q 01 +0.006 
u - ul -0.004 

0.01 ±0.01 






0.03 ±0.03 


- 7T°7T° 




< 0.01 


n + X -> n°X' 




04+ 04 


Total 


12 +(J - US 

u - iz -0.04 



TABLE II. Summary of the systematic error on the 
SES. We refer to the error coming from the statistics of 
Kl — > e + e~7 events as "statistical error", and the remain- 
ing errors as "systematic error". 



Source of Error 


Size(%) 


Statistical error 


0.79 


B(K L -> e + e"7) 


5.5 


B(n Q -> e + e- 7 ) 


2.7 


MC statistics 


0.81 


Drift Chamber inefficiency 


1.89 


Energy measurement 


0.80 


Momentum measurement 


0.07 


TRD efficiency 


1.88 


Total of systematic error 


6.78 
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